Obscured star formation in the central region of the dwarf galaxy 

NGC 5253 1 

Almudena Alonso-Herrero 2 
Steward Observatory, The University of Arizona, 933 N. Cherry Ave., Tucson, AZ 85721 

and 

Departamento de Astrofisica Molecular e Infrarroja, IEM, Consejo Superior de 
Investigaciones Cientificas, Serrano 113b, 28006 Madrid, Spain 

Toshinobu Takagi 

Centre for Astrophysics and Planetary Science, University of Kent, Canterbury, Kent, CT2 

7NR, UK 

Andrew J. Baker 

Max-Planck-Institut fur extraterrestrische Physik, Postfach 1312, D-85741 Garching, 

Germany 

George H. Rieke, Marcia J. Rieke 

Steward Observatory, University of Arizona, 933 N. Cherry Ave, Tucson, AZ85721 

Masatoshi Imanishi 2 

National Astronomical Observatory 2-21-1 Osawa, Mitaka, Tokyo 181-8588, Japan 

and 

Nick Z. Scoville 
California Institute of Technology, Pasadena, CA 91125 

ABSTRACT 



2 Visiting Astronomer at the Infrared Telescope Facility, which is operated by the University of Hawaii 
under Cooperative Agreement no. NCC 5-538 with the National Aeronautics and Space Administration, 
Office of Space Science, Planetary Astronomy Program. 



- 2 - 



We present JZST/NICMOS observations (1.1 — 2.2 /xm) and 1.9 -4.1 /wn spec- 
troscopy of the central region of the dwarf galaxy NGC 5253. The i/ST/NICMOS 
observations reveal the presence of a nuclear double star cluster separated by 
~ 0.3 — 0.4" or 6 — 8pc (for a distance d = 4.1 Mpc). The double star cluster, 
also a bright double source of Pact emission, appears to be coincident with the 
double radio nebula detected at 1.3 cm. The eastern near-infrared star cluster 
(CI) is identified with the youngest optical cluster, whereas the western star 
cluster (C2), although it is almost completely obscured in the optical, becomes 
the brightest star cluster in the central region of NGC 5253 at wavelengths 
longer than 2 /iin. Both clusters are extremely young with ages of approxi- 
mately 3.5 million years old. C2 is more massive than CI by a factor of 6 to 
20 (M C2 = 7.7 x 10 5 - 2.6 x 1O 6 M , for a Salpeter IMF in the mass range 
0.1 — 100 M ). Analysis of the circumnuclear spectrum excluding CI and C2, 
as well as of a number of other near-infrared selected clusters with a range of 
(young) ages, suggests that the star formation was triggered across the central 
regions of the galaxy. We have also modelled the nuclear UV to mid-infrared 
spectral energy distribution (SED) of NGC 5253 and found that the infrared 
part is well modelled with a highly obscured (Ay = 17 mag) young starburst 
with a stellar mass consistent with our photometric estimates for CI and C2. 
The SED model predicts a moderately bright polycyclic aromatic hydrocarbon 
(PAH) feature at 3.3 /xm that is not detected in our nuclear L-band spectrum. 
NGC5253's low metallicity and a top-heavy IMF likely combine to suppress the 
3.3 /im PAH emission that is commonly seen in more massive starburst systems. 

Subject headings: galaxies: individual: NGC5253 — galaxies: dwarf — galaxies: 
nuclei — galaxies: starburst — galaxies: star clusters 



1. Introduction 

It has long been recognized that the starburst occurring in the central region of NGC 5253 
must be among the youngest in the local Universe, and that this star forming process may be 
the first to have occurred in this galaxy in the last ~ 10 8 yr (van den Bergh 1980; Moorwood 



1 Based on observations with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope 
Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc. under 
NASA contract No. NAS5-26555. 



- 3- 



& Glass 1982; Rieke, Lebofsky, & Walker 1988). NGC 5253 appears to have been a metal 
poor dwarf elliptical before the current burst of star formation, and after the current burst 
of star formation it may evolve back to a dwarf elliptical galaxy with a bright core such as 
those in the Virgo and Fornax clusters (Caldwell & Phillips 1989). 

One of the longstanding questions is the triggering mechanism for the central starburst 
of NGC 5253. It has been suggested that the nuclear burst of star formation was triggered as 
the result of an interaction with the neighboring galaxy M83 (van den Bergh 1980; Caldwell 
& Phillips 1989). Recently Meier, Turner, & Beck (2002) have found that only the giant 
molecular cloud in the dust lane intersecting the galaxy, among the five detected in CO(2-l), 
is sufficiently nearby to be associated with the central starburst. More interestingly, these 
authors have shown that gas is falling into the nucleus of this galaxy, suggesting that this 
may be the mechanism for fuelling the starburst. 

Evidence for an extremely young and powerful burst of star formation in the central 
region of NGC 5253 comes from different observational clues. The nuclear radio and mil- 
limeter continuum emission are almost entirely due to thermal emission from Hll regions 
(Beck et al. 1996; Mohan, Anantharamaiah, & Goss 2001; Meier et al. 2002). Both the 
detection of Wolf-Rayet optical features (Walsh & Roy 1987; Conti 1991) and the high stel- 
lar temperature (T > 35, 000 — 40, 000 K), inferred from the near-infrared Hel(1.7yum)/Brl0 
line ratio (Vanzi & Rieke 1997) and mid-infrared emission lines (Beck et al. 1996; Crowther 
et al. 1999; Thornley et al. 2000), can only be explained with the presence of massive stars 
(IMF extending to at least M ~ 40 - 60 M , see Rigby & Rieke 2004). Finally Cowan 
& Branch (1982) searched for two historical supernovae in this galaxy using the VLA, and 
did not detect them, or any other radio supernovae, indicating that the burst must be very 
young. The star formation appears to be occurring in a relatively small region of the galaxy 
on the observed Ha emission (Calzetti et al. 1997) and the compact mid-infrared emission 
of this galaxy (Rieke 1976). Thus, NGC 5253 offers us an ideal opportunity to study the 
early stages of extragalactic star formation in the relatively confined nuclear region of this 
galaxy. 

Throughout the paper we assume a distance to NGC 5253 of d — 4.1 Mpc (Sandage et 
al. 1994). 

2. Observations 



2.1. iZST/NICMOS Near-Infrared Imaging 

#5T/NICMOS observations of NGC 5253 were obtained in January 1998 using the NIC2 
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camera, with a pixel size of 0.076" pixel" 1 . The filters used were F110W, F160W and F222M, 
roughly equivalent to ground-based filters J, H and K respectively. In addition, narrow band 
observations through the F187N and F190N filters, covering the Paa (A rcst = 1.87/xm) + 
continuum and 1.90 /im continuum, respectively, were obtained. 

The raw images were reduced with routines from the package NicRed (McLeod 1997). 
The main steps in the data reduction involve subtraction of the first readout, dark current 
subtraction on a readout-by-readout basis, correction for linearity and cosmic ray rejection 
(using FULLFIT), and flat fielding. The flux calibration was performed using the conversion 
factors derived from measurements of the standard star P330-E during the Servicing Mission 
Observatory Verification (SMOV) program (Marcia Rieke, private communication 1999). 

NICMOS only provides continuum bands to the red of the emission lines. The extinction 
in the nuclear region of NGC 5253 is known to be very high (even at infrared wavelengths; 
e.g., Turner et al. 2003 and references therein), and a straight subtraction of the longer 
wavelength continuum image may result in an overcorrection of the continuum at the wave- 
length of the emission line. We fit the continuum between 1.60 /im and 1.90 /im (using the 
NIC2 F160W and NIC2 F190N line-free images) on a pixel-by-pixel basis, and interpolate to 
estimate the continuum at 1.87/zm. This interpolated continuum is then subtracted from the 
line+continuum image to produce the continuum subtracted Pact image. The fully-reduced 
images were rotated to the usual orientation with north up, east to the left. They are shown 
in Fig. 1. In addition, we constructed infrared color maps equivalent to ground-based J — H 
(F110W - F160W) and H — K (F160W - F222M) color maps. 

2.2. IRTF/Spex Near-Infrared K and L-band Spectroscopy 

We have obtained simultaneous intermediate resolution K (A = 1.9 — 2.4 /im) and L 
band (A = 3 — 4.1 /im) spectroscopy of the central region of NGC 5253 using Spex (Rayner 
et al. 2003) at the 3 m NASA IRTF telescope on Mauna Kea in March 2003. We used 
the 0.15" pixel _1 plate scale and a slit width of 1.6" which provided a spectral resolution 
of R ~ 600 in the i^-band. The slit was oriented approximately North-South to cover the 
nucleus (C1+C2, Section 3.1, Fig. 2) of NGC 5253 plus the bright star cluster C3+C3' to 
the southwest of the nucleus (see Table 1 and Fig. 2). The C4+C5 cluster (see Fig. 2) was 
sufficiently bright and pointlike at the guider resolution that we could use it for the tip-tilt 
capability, thus improving the spatial sampling of the data. The total on-source integration 
time was 30 minutes. The seeing conditions as measured from standard star observations 
were approximately 0.6" (FWHM) in the i^-band. We also obtained KL spectroscopy of 
two comparison galaxies Henize 2-10 and NGC 2903 in a similar fashion. Observations of 
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solar-type stars were interspersed with the galaxy observations selected to be at similar air 
masses, and used to perform the telluric correction. 

We used the data reduction package SpexTool (Cushing, Vacca, & Rayner 2004) 
especially developed to reduce Spex data. Briefly, we created calibration frames (flatfields 
and arcs), median combined the galaxy observations, flatfielded the data, extracted ID 
spectra, and performed the wavelength calibration using sky lines and arcs. The 1-D spectra 
of the nuclear region of NGC 5253 were extracted with apertures: 1.5" x 1.6" and 3" x 1.6". 
We also extracted a i-T-band spectrum of the C3+C3' cluster through a 1.5" x 1.6" aperture. 

Standard stars at a similar airmass observed before and after the galaxy were reduced 
in a similar fashion and used to perform the correction of the galaxy spectra for atmospheric 
transmission. The flux calibration of the galaxy spectra was obtained from the observations 
of the standard stars. We compared the Spex .fT-band continuum level with the flux measured 
in the NIC2 F222M images through the same extraction apertures and found that they agree 
to within 10%. The fully reduced Spex spectra corrected for the standard star blackbody 
temperature, flux-calibrated and shifted to rest frame wavelengths are shown in Fig. 3. Line 
fluxes and equivalent widths (EW) for the main emission lines are given in Table 2. 

2.3. HST/WFPC2 UV and Optical Imaging 

Fully reduced HST/WFPC2 images of NGC 5253 through the F255W, F547M and 
F814W filters (Calzetti et al. 1997) were retrieved from the HST archive. These observations 
placed the central part of the galaxy on the WP3 chip providing a plate scale of ~ 0.1" 
pixel -1 . In addition we retrieved images through the F555W and F814W filters where 
the galaxy was centered on the WP1 chip providing higher spatial resolution. Because 
the [Olll]A5007A emission line in the central region of NGC 5253 is bright (e.g., Gonzalez- 
Riestra, Rego, & Zamorano 1987), we expect it to have a large contribution to the photometry 
through the F555W filter. We thus use the F555W image for comparing the optical and near- 
infrared morphologies, but not for star cluster photometry. 

2.4. ISO/ISOCAM Mid-Infrared Imaging 

We obtained fully reduced ISO/ISOCAM images of NGC 5253 observed through five 
filters: LW1 (4.5 H. L W2 (6.7 /im), LW6 (7.7/xm), LW10 (12 /W), and LW9 (14.9 /xm). 
The pixel size is 1.5" pixel -1 , which provides a field of view of 48" x 48". 
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Fig. 1.— HST/NICMOS images of the central 19" x 19" of NGC 5253. Orientation is North 
up, East to the left. The continuum images are the J-band image (NIC2 F110W filter) 
on the upper left panel, the if-band image (NIC2 F160W) on the upper right panel, and 
the if-band image (NIC2 F222M) on the lower right panel. The continuum-subtracted Paa 
(NIC2 F187N - F190N) line image is on the lower left panel. The double star cluster C1+C2 
(H1+H2) is the bright source on the upper left of each image (see also Fig. 2). 
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Fig. 2.— Close-ups of the central 11.9" x 11.9" of NGC 5253. Orientation is North up, East 
to the left. Upper left panel: Contours of the HST/WFPC2 F814W (WP1 chip) image. The 
pixels have been rebinned to match the pixel size of the NIC2 images. Upper right panel: 
J — H color map, where dark means redder J — H colors. Lower left panel: Contours of 
the if -band emission (NIC2 F160W). We mark the positions of the five brightest if-band 
clusters. Lower right panel: Contours of the continuum-subtracted Paa (NIC2 F187N - 
F190N) emission. 
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Fig. 3. — Upper panel: Spex if-band spectra of the nuclear region of NGC 5253 extracted through two 
different apertures: 1.5" x 1.6" and 3" x 1.6". In addition we show spectra of a circumnuclear region (that 
is, excluding CI and C2, see Section 5.1) and the bright cluster C3+C3' to the southwest of the nucleus (see 
Fig. 2 and Table 1 for location with respect to the nucleus of NGC 5253). On top of the plot we show the 
positions of emission lines positively identified, as well as tentative identifications of H2 lines at 1.967 /jm and 
2.1218 /im, and Hei lines at 2.1128 /im and 2.1137 /im (see Lumsden, Puxley, & Doherty 1994). We also mark 
the positions of the bandhead of the CO features (Kleinmann & Hall 1986). For display purposes the if -band 
spectra have been shifted along the vertical direction by the offsets indicated in the figure. Lower panel: 
Spex L-band spectra of the central region of NGC 5253 extracted through two different apertures: 1.5" x 1.6" 
and 3" x 1.6", and the circumnuclear region. The spectra have been shifted to rest-frame wavelengths. We 
also mark the positions of bright emission lines. 



Table 1. Clusters in the central region of NGC 5253 with M H < — 11 mag. 



Cluster 


X 


Y 


-MF160W 


A4"F222M 


m-FllOW — "JF160W m-Fl60W — m F222M 


logEIU(Paa) 


Possible ID 




(") 


(") 










(A) 




Central Double Cluster 


C2 


0.0 


0.0 


-11.68 ± 0.13 


-14.16 ± 0.06 


1.6 ± 0.1 


2.5 ± 0.1 


3.48 ± 0.23 


Dominant 1.3cm source 


CI 


0.3E 


0.1S 


-11.57 ±0.14 


-12.21 ±0.16 


0.3 ±0.2 


0.6 ±0.2 


3.41 ±0.29 


NGC5253-5, UV-12, Harrisl 


















Secondary 1.3cm source 












Other Brij 


;ht Clusters 






C3* 


0.1E 


7.0S 


-12.74 ±0.16 


-12.38 ±0.13 


0.9 ±0.2 


-0.4 ±0.2 




NGC5253-1, UV-2, Harris2 


C3'* 


0.3E 


7.2S 


-11.43 ± 0.11 


-12.72 ±0.10 


1.0 ±0.1 


1.3 ±0.2 






C3+C3' 


0.1E 


7.1S 


-13.48 ± 0.05 


-13.85 ±0.05 


1.0 ±0.1 


0.4 ±0.1 


1.26 ±0.20 


(0.3WJ.0S) 


C4 


5.7W 


4.9S 


-12.67 ±0.07 


-12.83 ±0.06 


0.9 ±0.1 


0.2 ±0.1 




NGC5253-3, UV-4, Harris3 


C5 


5.1W 


4.5S 


-11.85 ± 0.15 


-11.53 ±0.09 


0.8 ±0.2 


-0.3 ±0.2 


0.88 ±0.20 


NGC5253-2, Harris5 


C6 


0.4E 


2.5N 


-11.43 ± 0.10 


-11.61 ±0.06 


1.0 ±0.1 


0.2 ±0.1 




Harris32 


C7 


1.6W 


13.9S 


-11.32 ± 0.04 


-11.65 ±0.03 


1.0 ±0.1 


0.3 ±0.1 




Harris27 


C8 


7.9W 


9.2S 


-11.19 ± 0.08 


-11.41 ±0.06 


0.7 ±0.1 


0.2 ±0.1 




NGC5253-6?, Harris9 


C9 


1.3W 


4.8S 


-11.10 ± 0.03 


-11.36 ±0.02 


1.0 ±0.1 


0.3 ±0.1 


1.47 ±0.09 




CIO 


4.0W 


12.8S 


-11.07 ±0.05 


-11.41 ±0.04 


0.9 ±0.1 


0.3 ±0.1 






Cll 


6.1W 


10.3S 


-11.06 ± 0.06 


-11.39 ±0.05 


1.0 ±0.1 


0.3 ±0.1 






C12 


7.8W 


9.0S 


-11.04 ±0.08 


-11.25 ±0.05 


0.9 ±0.1 


0.2 ±0.1 


1.68 ±0.18 


NGC5253-6? 



Note. — X and Y are the measured offsets in arcscconds relative to the position of C2. The magnitudes, colors, and EWs are 
measured through a 0.76"-diameter aperture (except for C3±C3', see below), and are not corrected for extinction. The column 
'ID' indicates possible correspondences with sources detected at other wavelengths. "NGC5253-1, NGC5253-2, NGC5253-3, 
NGC5253-5, and NGC5253-6" arc optical star clusters from Calzetti et al. 1997. "UV-2, UV-4 and UV-12" are UV clusters 
identified by Meurer et al. 1995. Also listed are correspondences with clusters reported by Harris ct al. (2004). The 1.3cm 
radio sources are from Turner et al. 2000. The near-infrared hot spots in Davics ct al. 1998 are named after the offsets quoted 
in their paper relative to the nucleus, where (0.3W,7.0S) corresponds to C3±C3' and (5.8W,4.6S) to C4+C5. 
"The C3 and C3' clusters are severely blended and thus the colors and magnitudes are strongly affected by deconvolution errors. 
Wc give magnitudes and colors for C3±C3' through a 1.5"-diamcter aperture. 
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Table 2. Line fluxes: //5T/NICMOS Paa measurements, and Spex K- and L-band 
spectroscopy of the nuclear region of NGC 5253. 



Line 


^rest 


Aperture 


Flux 


EW 








/ 9 1 \ 

(erg cm s J 


(A) 


Paa 


1.877 


1.5" x 1.6" 


4.75 x 10~ 13 


~ 2000 






3" x 1.6" 


5.96 x 10~ 13 


~ 2200 


Bi5 


1.945 


1.5" x 1.6" 


3.40 x 10~ 14 


177 






3" x 1.6" 


4.60 x 10~ 14 


146 


Hel 


2.058 


1.5" x 1.6" 


2.48 x 10- 14 


113 






3" x 1.6" 


3.75 x 10~ 14 


116 


Br7 


2.166 


1.5" x 1.6" 


5.72 x 10~ 14 


255 






3" x 1.6" 


8.25 x 10~ 14 


250 






circumnuclear 


2.67 x 10~ 14 


216 






cluster C3+C3' 




< 5 


Pf 7 


3.741 


1.5" x 1.6" 


4.43 x 10~ 14 


32 






3" x 1.6" 


5.90 x 10- 14 


37 


Bra 


4.051 


1.5" x 1.6" 


3.62 x 10~ 13 


205 






3" x 1.6" 


4.64 x lO" 13 


235 



Note. - Paa measurements are from the narrow-band iJST/NICMOS imaging data. 
Other measurements are from the IRTF Spex spectroscopy. 
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Fig. 4. — Close-ups of the central 3" x 3" of NGC 5253, showing the two nuclear com- 
pact Hll regions (Hl+H2)/star clusters (C1+C2). The thin lines are the contours of the 
continuum-subtracted Paa emission (left panel, in red in the electronic edition), and the 
1.6 jum continuum emission (right panel, in blue in the electronic edition). In both panels 
the thick contours are the 1.3cm VLA radio map from Turner et al. (2000) showing the 
location of the two radio sources. Orientation is North up, East to the left. 
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Fig. 5.— False color image produced with the HST/WFPC2 F555W, iJST/NICMOS F110W, 
and #ST/NICMOS F160W images. The FOV is similar to that in Fig. 2. Note the red color 
of the C2 cluster. 



-13- 



3. Photometry 

3.1. iJ-band selected star clusters 

We choose to identify the star clusters in the central region of NGC 5253 using the 
if-band (NIC2 F160W filter) data because this wavelength represents the best compromise 
between spatial resolution and reduced effects of extinction. We used IRAF daophot rou- 
tines and detected a total of 269 clusters. The positions of the clusters identified in the 
if-band were then used as input for daophot at the other four wavelengths, 1.1, 1.87, 1.9 
and 2.2 /im. Since most of the clusters appear to be slightly resolved, we constructed at each 
wavelength a point spread function (PSF) from bright sources in the images to model the 
photometry instead of using theoretical PSFs. These PSFs are used by daophot to model 
the cluster photometry simultaneously for all the sources detected in the field of view. The 
cluster photometry was obtained through a 10 pixel=0.76" (~ 15 pc) diameter aperture. The 
underlying galaxy emission was measured from annuli surrounding the source. In addition 
to the continuum photometry, we compute for each cluster the EW of Paa. 

Table 1 lists the photometry for the brightest star clusters selected in the if -band with 
absolute magnitudes M F160W < — 11. The quoted errors are the output of DAOPHOT. Since 
most of the clusters, and in particular the bright ones, appear to be slightly resolved, we did 
not apply an aperture correction to the magnitudes measured through the 0.76"-diameter. 
The aperture corrections (for point sources) are —0.18, —0.19 and —0.28 for the F110W, 
F160W and F222M magnitudes, respectively, and +0.09 mag for the mpieow — m F222M colors, 
based on photometry of point sources generated with Tiny Tim (Krist et al. 1998). 

3.2. H II regions 

We identified and performed photometry of Hll regions in the continuum subtracted 
Paa image of the central 19" x 19" region of NGC 5253. We used the software region, 
kindly provided by C. H. Heller (see Pleuss, Heller, & Fricke 2000, and references therein for 
a detailed description). 

region is a semiautomated method that locates and computes statistics of H II regions 
in an image, based on contouring taking into account the local background. We set the 
limit for the size of an H II region to 9 contiguous pixels, which at the distance of NGC 5253 
corresponds to minimum linear sizes (diameters) of approximately 5pc. Briefly, the Hll 
selection and photometry procedure is as follows. First the program selects local maxima; 
then we set a brightness limit for the pixel to be included as part of an Hll region that 
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the pixel must have an intensity above the local background of at least 3 times the rms 
noise of the local background (see Rand 1992 for more details on the criteria employed). 
After identifying the Hn region, the software measures its position, size (area), and Pact 
flux by subtracting the closest local background (selected interactively by the user) from the 
observed flux. 

In Table 3 we present the photometry of H II regions with measured Pact luminosities 
(not corrected for extinction) log (Lp aQ ,/ergs _1 ) > 38. The diameters listed in this table 
are 'equivalent diameters' computed assuming that the area covered by the Hn region is 
circular. 



3.3. Nuclear Spectral Energy Distribution 

In deriving the nuclear spectral energy distribution (SED) we can only zoom in as 
tightly as the smallest of our extraction apertures (1.5" x 1.6") for the .fTL-spectroscopy. We 
thus matched the photometry at other wavelengths to this size, and paid special attention 
to subtracting the underlying galaxy emission to isolate as much as possible the nuclear 
emission. In what follows, we describe the procedure for obtaining the photometry and 
estimating the galaxy contribution. The nuclear photometry is summarized in Table 4. 

For the photometry of the HST/WFPC2 and NICMOS broad-band images we simulated 
the spectroscopy aperture on the images, and the local background was estimated from annuli 
around the nucleus taking care that they did not include faint star clusters. We estimated 
the continuum at 1.90 /im from the narrow-band F190N image in a similar fashion. 

The ISOCAM images of NGC 5253 show at all five wavelengths a bright point source 
superimposed on a faint extended component. The emission from the point source is unre- 
solved, as determined from the comparison of the measured FWHM of the nuclear source 
with the sizes of theoretical PSFs (with sizes typically of 2 — 4.5" FWHM depending on the 
wavelength, Okumura 1998). We assume that the unresolved flux at all five wavelengths 
arises from the nuclear region. The local background was estimated from areas near the 
edges of the images. 

We also make use of mid-infrared data from the literature. The mid-infrared narrow- 
band and broad-band fluxes of Frogel, Elias, & Phillips (1982) were observed through an 
8"-diameter aperture. This photometry is in good agreement with the 8 — 13 /im spectropho- 
tometry through a 5.4"-diameter aperture of Aitken et al. (1982). Other ground-based 
measurements from the literature include the 11.7 and 18. 7 //m unresolved fluxes from Gor- 
jian et al. (2001), and the 10 and 21//m fluxes through a 5.4"-diameter aperture of Rieke 
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(1976) and Lebofsky & Rieke (1979). From comparison of the ground-based measurements, 
the only discrepant point for the nuclear photometry is the 8" -diameter flux at 20 /im from 
Frogel et al. (1982). The ground-based small aperture mid-infrared fluxes are consistent 
with the ISOCAM unresolved fluxes, giving us confidence that most of the small aperture 
mid-infrared emission is from the nucleus. 2 The IRAS fluxes (Roche et al. 1991) and the 
SCUBA (James et al. 2002) data point will be taken as upper limits for the modeling of the 
nuclear SED (Section 7) as they encompass the entire galaxy. 

4. Morphology of the central super star cluster/Hn Region 

The most striking feature revealed by the HST/NICMOS images is the presence of a 
double star cluster (see also Turner & Beck 2004) in the nucleus of NGC 5253, which we 
term CI for the eastern cluster and C2 for the western cluster (see Table 1). Fig. 2 shows 
close-ups of the central 11.9" x 11.9" of NGC 5253 in the optical (WFPC2 F814W) and at 
near-infrared wavelengths. We also indicate the locations of the five brightest clusters in the 
-ff-band (see also Table 1). In Fig. 4 the close-ups show the central 3" x 3". 

The double star cluster is most conspicuous at 1.6 /im (see Fig. 4), where both clusters 
show similar apparent brightness (see Table 1). The distance between the two star clusters 
as measured at 1.6 /tm is 0.3 — 0.4" or 6 — 8pc for the assumed distance. Turner et al. 
(2003) have recently reported that the brightest source at 2.2 /im is offset by ~ 0.3 ± 0.1" 
to the northwest of the youngest optical source and peak of the Ha emission — referred to 
as NGC 5253-5 and UV-12 in Calzetti et al. (1997) and Meurer et al. (1995), respectively. 
Indeed the NICMOS F222M image shows that the western cluster C2 dominates the emission 
at 2.2 /im. We thus identify the eastern star cluster CI with the peak of the Ha emission and 
youngest optical star cluster (NGC 5253-5 in Calzetti et al. 1997 notation). The western 
star cluster C2 becomes the brightest source at 2.2 /im (see Table 1, and Fig. 2), and is so 
reddened that it is hardly detected in the optical (see the WFPC2 F814W image shown in 
Fig. 2). This fact is also illustrated in the false color image of the central region of NGC 5253 
shown in Fig. 5. 

Bright Paa emission is detected from both star clusters, and the two Hn regions are 
named H2 (at the location of C2) and HI (at the location of CI), see Table 3. The Ha 



2 We note that Rieke (1976) and Lebofsky & Rieke (1979) already remarked that NGC 5253 has a compact 
10 /xm source with a much higher surface brightness than other star forming galaxies in their sample. Gorjian 
et al. (2001) used Keck and measured a size for the central mid-infrared source of 0.62", supporting the fact 
that the mid-infrared flux originates in a very small nuclear region. 
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luminosity derived from the Paa luminosity (not corrected for extinction) using the theoret- 
ical line ratios in Hummer & Storey 1987) of H2 indicates that this source has the typical 
brightness of giant Hn regions, e.g., 30 Doradus. 3 

Turner, Beck, & Ho (2000) obtained 1.3 cm and 2 cm VLA images of NGC 5253 at 
subarcsecond resolution and detected at both wavelengths a bright supernebula. Their 1.3 cm 
map also displays a secondary peak within the central 0.5". Fig. 4 compares the radio and 
near-infrared morphologies (Paa emission and 1.6 fim) of the central 3" x 3". Assuming that 
the brightest radio source and the C2 star cluster (also the H2 Hn region) are coincident, 
then the secondary radio peak appears to be at the location of secondary Hn region (HI) 
and at the approximate location of the star cluster CI. The central double cluster/Hn region 
appears also to be coincident with the bright mid-infrared source detected by Gorjian, Turner, 
& Beck (2001). It is likely that the bright near-infrared cluster C2 is actually responsible for 
most of the mid-infrared emission (see Turner & Beck 2004). 



3 Note that the near-infrared magnitudes and EWs of CI and C2 do not correspond to exactly the 
same spatial extent as the double H n region because of the two (unavoidably) different methods for doing 
photometry of star clusters (approximately point sources) and H II regions (extended sources) . 
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Table 3. Hn regions in the central region of NGC 5253 with log (L PaQ /ergs x ) > 38. 



HII 


X 


Y 


Diameter logL(Paa) 


ID 


Region 


(") 


(") 


(") (erg s- 1 ) 










Central Double Hn Re 


gion 


H2 








2.0 38.91 


C2, Dominant 1.3cm source 


HI 


0.3E 


0.1S 


1.1 38.42 


CI, Secondary 1.3cm source 










NGC5253-5, UV-12 


Other H II Regions 


H3 


0.8W 


0.1S 


1.2 38.11 




H4 


LIE 


1.7S 


1.5 38.07 




H5 


0.8W 


1.1N 


1.1 38.00 




H6 


0.7W 


2.6S 


1.3 38.00 


NGC5253-4, UV-1, Harris4 


Note. - 


- X and Y are 


the measured offsets in 


arcseconds relative to the position of H2. 



The Pact luminosities are not corrected for extinction. The diameters listed in this table are 
'equivalent diameters' computed assuming that the covered area by the H II region is circular. 
The Pact luminosities correspond to these areas, and are not corrected for extinction. The 
column 'ID' indicates possible correspondences with sources detected at other wavelengths. 
"NGC5253-5/UV-12" and "NGC5253-4/UV-1" are young star clusters detected in both the 
optical and UV, Calzetti et al. 1997 and Meurer et al. 1995, respectively. "Harris4" is an 
optical cluster reported by Harris et al. (2004). The 1.3cm radio sources are from Turner et 
al. 2000. 



-18- 



Table 4. Photometry of nuclear region of NGC 5253 



A 


J? 

Jv 


Aperture 


Reference 


Uncertainty 


fj,m 


mjy 








0.255 


0.20 


1.5"xl.6" 


WFPC2 F255W 


7% 


0.547 


0.40 


1.5"xl.6" 


WFPC2 F547M 


3% 


0.814 


0.84 


1.5"xl.6" 


WFPC2 F814W 


5% 


1.1 


1.20 


1.5"xl.6" 


NICMOS NIC2 F110W 


7% 


1.6 


1.20 


1.5"xl.6" 


NICMOS NIC2 F160W 


6% 


1 Q 

i.y 


l.Oo 


1.0 XI. 


lNl^lvlWo INIOZ r lyuiN 


/o 


2.2 


3.55 


1.5"xl.6" 


NICMOS NIC2 F222M 


4% 


z.z 


a on 


1.0 XI. 


Spcx spectroscopy 


lu/o 


3.0 


19.0 


1.5"xl.6" 


Spcx spectroscopy 


10% 


3.5 


48.4 


1.5"xl.6" 


Spcx spectroscopy 


10% 


4.0 


92.5 


1.5"xl.6" 


Spcx spectroscopy 


10% 


7.8 


867. 


circular-8" 


Frogel et al. 1982 


27% 


O.D 


SI 7 
ol / . 


circular-8" 


rrogei er ai. iyoz 


lo/o 


9.6 


1067. 


circular-8" 


Frogel et al. 1982 


14% 


10.4 


1541. 


circular-8" 


Frogel et al. 1982 


11% 


11.4 


1574. 


circular-8" 


Frogel et al. 1982 


12% 


12.4 


2041. 


circular-8" 


Frogel et al. 1982 


22% 


10.0 


1298. 


circular-8" 


Frogel et al. 1982 


6% 


20.0 


6098. 


circular-8" 


Frogel et al. 1982 


15% 


10.0 


1500. 


circular-5.4" 


Lebofsky & Rickc 1979 




21.0 


2800. 


circular-5.4" 


Lcbofsky & Rickc 1979 


25% 


11.7 


2200. 


unresolved 


Gorjian et al. 2001 


10% 


18.7 


2900. 


unresolved 


Gorjian et al. 2001 


10% 


4.5 


143. 


point-source 


ISOCAM 


4% 


6.7 


520. 


point-source 


ISOCAM 


9% 


7.7 


829. 


point-source 


ISOCAM 


7% 


12.0 


2060. 


~point-source 


ISOCAM 


4% 


14.9 


4100. 


~point-source 


ISOCAM 


4% 


12.0* 


2730. 


entire galaxy 


IRAS-Roche ct al. 1991 




25.0* 


12630. 


entire galaxy 


IRAS-Roche et al. 1991 




60.0* 


30630. 


entire galaxy 


IRAS-Rochc ct al. 1991 




100.0* 


29490. 


entire galaxy 


IRAS-Roche et al. 1991 




850.0* 


192. 


41" 


SCUBA-Jamcs ct al. 2002 


12% 



Note. — * Upper limits to the nuclear emission. 
The errors account for the photometric calibration uncertainty, and when pos- 
sible, the uncertainty from estimating the contribution from the underlying 
galaxy (see Section 5.1). The absolute photometric calibration uncertainties 
for the ISOCAM data are from the ISOCAM handbook (Blommaert ct al. 
2001). We also make use of the 8 — 13 fim spectrophotometry of Aitkcn ct al. 
(1982). 
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Table 5. Extinction to the gas for the central region of NGC 5253 for a foreground dust 

screen model. 



Aperture 


Lines 


Ay 


Ref. 


Infrared estimates 


0.6" x 3" 


Bra/Br7 


18 ± 5mag 


2 


1.5" x 1.6" 


Bra/Br7 


13 ± 2mag 


1 




Paa/Br7 


13 ± 5mag 


1 




Br5/Br7 


6 ± 5mag 


1 


10.3" x 20.7" 


Bra /Br 7 


8 ± 3 mag* 


3 


Circular 5.4" 


9.7/im Silicate feature 


lOmag 


3,4 


Optical estimates 


Central 23" x 23" 


Ra/R(3 


1-3 mag** 


5 



Note. — References: 1. This work. 2. Turner et al. 2003. 3. Kawara, Nishida, & Phillips 
1989. 4. Aitken et al. 1982. 5. Walsh & Roy 1989. 

* Ay obtained from the Bra/Br7 line ratio reported in Kawara et al. 1989 and the extinction 
of Rieke & Lebofsky 1985 

** Observed range of Ay in the central region of NGC 5253. 
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5. The Nuclear Stellar Population 
5.1. Age 

The spectroscopy through the 1.6" x 1.5" aperture can be used to derive the properties of 
the stellar population of the central cluster C1+C2, as long as it is not heavily contaminated 
by other faint clusters/Hn regions. We can compare the Pact luminosities of the Hll regions 
HI and H2 (Section 3.2) believed to represent well the luminosities of the two nebulae, with 
that measured simulating the Spex 1.6" x 1.5" aperture on the continuum subtracted Pact 
image. We find that these two measurements are within ~ 15% of each other, and thus we 
believe that the spectroscopic data are mainly sampling the double cluster/Hn region. 

The EWs of hydrogen recombination lines are reliable age estimators of young (that 
is, the lines are in emission) stellar populations for the case of instantaneous star formation 
(e.g., Alonso-Herrero et al. 1996) since EWs measure the ratio between massive (young) 
stars and the total number of stars. Moreover, if the extinction to the gas is the same as 
the extinction to the stars (cf. Calzetti 2001) then the EWs are extinction-independent age 
indicators. The observed EW of Br7 (Table 2) yields an age for the nuclear stellar population 
of C1+C2 of approximately 3.4 ± 0.2 x 10 6 yr, similar to that estimated by Davies, Sugai, 
& Ward (1998) using Br7 narrow-band imaging. We have used Starburst99 (Leitherer et 
al. 1999) for a stellar population formed in an instantaneous burst with solar metallicity 
and a Salpeter IMF with lower and upper mass cutoffs of M\ ow = 1 M Q and M up = 100 M , 
respectively. This age estimate is not very dependent on the metallicity of the stars, but it 
is rather sensitive to the choice of IMF (see Leitherer et al. 1999). 

Using daophot on the HST/NICMOS images we have been able to separate both clus- 
ters/Hll regions and measure their near-infrared magnitudes and EW of Pact (see Table 1). 
Assuming that the age of young clusters is well represented by the observed EW of Pact, 
we find that CI and C2 are both young star clusters, with similar ages of approximately 
3.3 ± 1.0 x 10 6 yr, in good agreement with other estimates. For instance, Calzetti et al. 
(1997) obtained an age of 2.2 — 2.8 x 10 6 yr for the peak of the Ha emission (presumably our 
cluster CI) based on UV and optical photometric data, and the EW of Ha. Tremonti et al. 
(2001) using UV spectroscopy fitted the age of the brightest optical cluster to 2+^ x 10 6 yr. 
Ground-based and ISO mid-infrared spectroscopy of the central region of NGC 5253 suggests 
a stellar population 2 - 4 x 10 6 yr old (Beck et al. 1996 and Crowther et al. 1999). 

The non-detection of the if-band CO bands (see Fig. 3) sets an upper limit to the age of 
the population in the central double cluster of 7 x 10 6 yr. We also extracted a circumnuclear 
.fT-band spectrum by subtracting the 1.5" x 1.6" spectrum from the 1.5" x 3" one. Even in 
the region surrounding CI and C2 there is no evidence for the presence of CO bands. One 
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possibility for the non-detection of the CO bands even in the circumnuclear region would 
be the presence of dust emission peaking in the mid- infrared. This component would fill 
in the CO bands in the i^-band, and weaken them as observed in many Seyfert 1 and 2 
galaxies with a strong thermal infrared continuum (e.g., Ivanov et al. 2000). Indeed from 
Fig. 3 it is clear that in the nuclear region the 2 — 4 /itn continuum is very steep and may 
have a significant contribution from thermal emission (see discussion in Sections 5.2.2 and 
7.2). This is not however the case for the circumnuclear spectrum where the continuum slope 
resembles that of less reddened starbursts (see Fig. 3). Thus we believe that the lack of CO 
bands in the nuclear region of NGC 5253 is not fundamentally due to the presence of a red 
continuum. 

We note however that age estimates based on the CO bands may be dependent on the 
metal abundance ([Fe/H]) of the stars as there is a dependence of the strength of the near- 
infrared CO bands on the metallicity of the stars (Origlia et al. 1997; Frogel et al. 2001; 
Ivanov et al. 2004). The upper limit for the age of the stellar population based on the 
non-detection of the CO bands remains approximately constant (7 — 9 x 10 6 yr) even if we 
consider the models with the lowest metallicity in Starburst99 (Z = ^Z ). This together 
with the measured EW of Br7 of the circumnuclear region (Table 2) suggests that the stellar 
population within the entire central 3" x 1.6" (32 pc x 60 pc) region of NGC 5253 is very 
young (~ 3.5 x 10 6 yr), and it is not only confined to the double star cluster C1+C2. 

5.2. Extinction 

The interpretation of the properties of the central double cluster requires an accurate 
estimate of the extinction to the source. The extinction to the gas can be estimated from 
hydrogen recombination line ratios, so long as the comparison is based on measurements 
with the same beam sizes. To determine the extinction to the stars we can compare observed 
colors (or the SED) with the outputs of Starburst99 for the age range derived in the preceding 
section. The determination of the extinction to the gas and stars is described in the next 
two sections. 



5.2.1. Extinction to the gas 

There have been a number of estimates of the extinction to the gas in the central 
region of NGC 5253, based primarily on optical and near-infrared hydrogen recombination 
lines. The largest uncertainty in the extinction by far comes from the unknown distribution 
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of dust within the source. Meurer et al. (1995) studied the UV properties of a sample 
of starburst galaxies that included NGC 5253, and concluded that the dust geometry can 
be effectively described in terms of a foreground screen configuration near the starburst. 
However, their galaxies were UV-selected and may be biased toward low extinction and 
possibly an unrepresentative distribution of the dust. 

We use the observed emission line fluxes of Paa, BrS, Bvy and Bra and the extinction 
law of Rieke & Lebofsky (1985) to infer the extinction to the central 1.5" x 1.6". In Table 5 we 
show our results together with estimates from the literature using a foreground dust screen 
model. All the near-infrared estimates of the visual extinction to the gas are consistent 
to within the errors, whereas the optical estimates are significantly lower (see discussion in 
Calzetti et al. 1997). The differing values of the extinction as a function of wavelength are 
the usual indication that the extinction is either patchy or mixed with the emitting sources 
along the line of sight (or both). Indeed, Calzetti et al. (1997) have argued that a young 
star cluster will be deeply embedded in its natal molecular cloud and thus the stars, gas and 
dust are likely to be mixed, providing an efficient way of "hiding" the dust. 

Using a value of the extinction of Ay = 11 mag (which is consistent with all our esti- 
mates) and the observed fluxes through the 3" x 1.5" apertures we infer a number of ionizing 
photons iVLy = 4 x 10 52 s _1 . This is similar to, although slightly lower than, the radio esti- 
mates of Turner et al. (2000, 2003) and Turner & Beck (2004). The radio H92o; observations 
of Mohan et al. (2001) predict N Ly — 2 — 4 x 10 52 s~ 1 , but this line is only sensitive to a 
narrow range of nebular densities. 

5.2.2. Extinction to the stars and dust emission at near-infrared wavelengths 

A qualitative description of the extinction to the stars can be obtained from color maps 
(see the J — H color map of the central 12" x 12" in Fig. 2) and false color images (e.g., 
Fig. 5, and also Turner et al. 2003). The C2 cluster is heavily obscured, but it becomes the 
brightest source in this galaxy at 2.2 /xm. On the other hand, the region surrounding C2 and 
the region to the southeast of C2 appear blue with respect to other star clusters due to the 
youth of the ionizing stars. Also, there is evidence for the dust lane that crosses the galaxy 
approximately east-west. 

Evidently the high value of the central extinction derived in the previous section must be 
mostly associated with C2, as the infrared colors of CI appear to be quite blue (see Table 1) 
consistent with outputs from Starburst99 for a young starburst. A stellar population created 
in an instantaneous burst after ~3-3.4x 10 6 yr shows infrared colors: J — H = 0.12 — 0.04 
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and H — K = 0.35 — 0.19. The observed near-infrared colors of CI are consistent with 
Ay — 2 ± 2 mag. The J — H color of C2 is consistent with Ay ~ 14 ± 2 mag for the same age 
range as above. Using a model where the central cluster is hidden behind a thin dust layer 
and is completely embedded in a thick dust cloud Calzetti et al. (1997) found that values of 
the extinction Ay > 9 mag were necessary to account for the observed UV-optical colors of 
the central cluster. The dereddened H — K color of C2 is still very red. As we shall see in 
Section 7, this is the effect of hot dust emission with a significant contribution in the -fT-band 
and/or differential extinction. If we assume that there is no differential extinction between 
the H and the X-band, then approximately 70% of the i^-band emission in C2 could be 
produced by hot dust (with some small contribution from gas emission as well). 



5.3. Stellar mass of the double cluster 



Using the age and the extinction estimates, we can infer the stellar masses of the two 
clusters for a Salpeter IMF with mass cutoffs of M low = 1 M and M up = 100 M , using their 
absolute if -band magnitudes. For CI we estimate a stellar mass of M C i ~5x 10 4 M o (see 
also Fig. 6). If we continued the high mass slope all the way down to 0.1 M the total stellar 
mass would be Mqi — 1.3 x 1O 5 M . For the C2 cluster since the current dataset does not 
allow us to distinguish between hot dust emission in the JT-band and differential extinction, 
we can only obtain a range of stellar masses. If there is no differential extinction between the 
gas and the stars, and up to 70% of the if-band emission is due to hot dust, then the stellar 
mass of C2 would be M C2 > 3 x 10 5 M (or M C2 > 7.7 x 10 5 M for a Salpeter IMF extending 
down to 0.1 M ). The upper limit to this estimate would be if all the if-band emission were 
stellar in origin (implying that there is differencial extinction) and we corrected the observed 
H — K color to that corresponding to 3 — 3.4 Myr; in this case the stellar mass would be 
M C2 < 1 x 1O 6 M (or M C 2 < 2.6 x 10 6 M for an IMF down to 0.1 M Q ). In any case, C2 
dominates CI by mass for any plausible C2 age and extinction. Using the number of ionizing 
photons (corrected for extinction) derived in Section 5.2.1 we would require a mass for the 
two ionizing clusters of M C1+C2 ~lx 1O 6 M (or M C1+C2 = 2.6 x 1O 6 M for an IMF down 
to 0.1 M ), for the same star formation history. 
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6. ii-band selected star clusters in the central ~ 20" 

6.1. Bright clusters 

In Table 1 we give the photometry, colors and observed EW of Paa of bright clusters 
(M H < — 11 mag) selected in the if -band. In the last column of this table we give possible 
associations with star clusters detected at other wavelengths. The correspondences are solely- 
based on the measured relative offsets taking into account the fact that CI is easily identified 
in the UV and optical, and at near-infrared wavelengths A < 2 /im. 

We find near-infrared counterparts to five of the six star clusters optically selected by 
Calzetti et al. (1997), as well as some of the clusters reported by Harris et al. (2004), 
as shown in Table 1. The absolute Mpi6ow and Mp222M magnitudes (not corrected for 
extinction) are typical of near-infrared star clusters detected in other hot-spot galaxies and 
starbursts (e.g., Alonso-Herrero, Ryder, & Knapen 2001; Elmegreen et al. 2002; Maoz et 
al. 2001). Except for star cluster C2, the rest of the central clusters in NGC 5253 are not 
as bright as the near-infrared super star clusters identified in luminous and ultraluminous 
infrared galaxies (Scoville et al. 2000, Alonso-Herrero et al. 2002) and interacting galaxies 
(e.g., the Antennae, Gilbert et al. 2000; Mengel et al. 2002). 

At ground-based resolution the star cluster C3+C3' appears as the second brightest 
near-infrared cluster after C1+C2, and was identified by Davies, Sugai, & Ward (1998) 
as one of the bright three near-infrared hot spots (see Figs. 1 and 2) in this galaxy. The 
third brightest cluster identified in Davies et al. (1998) corresponds to the C4+C5 cluster 
in our notation (see Table 1). The if-band spectrum of C3+C3' is presented in Fig. 3. 
Note that the 1.5" x 1.6" extraction aperture possibly includes other fainter star clusters 
(see e.g., Fig. 2). The measured EW of Br7 sets a lower limit to the age of this cluster of 
approximately > 7 x 10 6 yr, using outputs of Starburst99 under the same assumptions used 
for the central cluster. This age is consistent with that derived from the observed EW of 
Pa« (7.5 ± 0.5 x 10 6 yr). The lack of 2.3 /im CO bands, on the other hand, sets an upper 
limit to the age of the cluster of approximately 7 — 9 million years, although this limit is 
dependent on the metallicity of the stars. Davies et al. (1998) using the measured EW of 
Br7 (= 14 — 16 A; note that they used a 2"-diameter aperture) from narrow-band imaging 
derived an age for this cluster of 10 million years, although using Starburst99 we would get 
~ 6.5 x 10 6 yr. Calzetti et al. (1997) estimated an age for C3 of 8 - 12 x 10 6 yr. 

Tremonti et al. (2001) studied in detail clusters C4 and C5 (clusters NGC5253-3 and 
NGC5253-2 in their notation) using HST UV spectroscopy, and derived ages of 3 Myr and 
8Myr, respectively. Their result is in good agreement with the age from the measured 
EW of Pact of C5 (see also Harris et al. 2004). However the young age Tremonti et al. 
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(2001) inferred for C4 appears to be inconsistent with the lack of bright Paa emission at the 
location of this star cluster. Harris et al. (2004) on the other hand find an age for this cluster 
of between 10 and 14 Myr consistent with our finding. For other near-infrared clusters in 
common with those optically selected by Harris et al. (2004) we find ages and stellar masses 
(see next section) consistent with their photometric and EW of Ha based estimates. 
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Young clusters (age<7xl0 yr) 




M F160W ( ma g) 

Fig. 6. — Main Panel: — Absolute Mpi6ow magnitude vs. equivalent width of Paa emission 
line for H-b&iad selected clusters (filled dots). We plot young star clusters with observed 
log EW(Paa)/A > 1.5 , that is, star clusters younger than approximately 7 x 10 6 yr. We mark 
the positions of the central double cluster CI and C2 (see Section 4.2) and Hn region H6 (see 
Section 6.2). The magnitudes have not been corrected for extinction. The arrow shows what 
would be the effect of correcting the observed absolute M F i6ow f° r 2 magnitudes of visual 
extinction assuming the same extinction to the stars and the gas. The lines represent the 
time evolution of star clusters with masses M = 10 3 M and M = 10 4 M for instantaneous 
star formation, a Salpeter IMF (Mi ow = 1M and M up = 100 M ), and solar metallicity 
using Starburst99 (Leitherer et al. 1999). The crosses on the model lines are drawn at 
intervals of 1 million years, with the youngest age on the upper part of the curves. Insert: 
- The empty histogram shows the distribution of the absolute Mpi6ow magnitudes of all 
the star clusters selected in the if-band, whereas the filled histogram shows the clusters with 
ages younger than 7 million years, that is, those shown in the main panel. 
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Fig. 7. — SED of the nuclear emission of NGC 5253. The small aperture measurements 
are from this work, Gorjian et al. (2001) and Lebofsky & Rieke (1979). As discussed in 
Section 3.3 the 20 jum flux from Frogel et al. (1982) is not used. For the nuclear SED 
the IRAS and Scuba data points are plotted as upper limits because these measurements 
comprise the entire galaxy. The best fit to the nuclear SED of NGC 5253 is the dashed line, 
whereas the dashed-dotted line shows the stellar contribution to this fit. For details on the 
modelling see Section 7. 
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6.2. Young versus old star clusters 

The only optically bright star cluster that is not bright in the near-IR is NGC5253-4, 
using the notation of Calzetti et al. (1997). This star cluster, however, is detected as a 
bright Hn region (H6, see Table 3) and is also the brightest UV source (Meurer et al. 1995). 
Based on the observed EW of Paa we derive a very young age for this cluster/Hll region 
(~ 3.5 x 10 6 yr), similar to that of CI and C2. Harris et al. (2004) derive an age for H6 
(their cluster # 4) of 5 x 10 6 yr using the EW of Ha. The youth of H6 explains why the 
ionizing cluster has not become infrared bright (that is, its absolute magnitude is fainter 
tha M F i6ow = — 11 mag). Starburst99 suggests that if the Hn region H6 was formed in an 
instantaneous burst the ionizing cluster will become near-infrared bright in the next three 
or four million years (see Fig. 6 for the location of H6 and future time evolution of its near- 
infrared magnitude). The age derived from the EW of Paa is consistent with that derived 
in Calzetti et al. (1997) using UV and optical magnitudes, together with the EW of Ha. 

As we have shown, the EW of Paa is a useful indicator of the age of young star clusters 4 . 
For all the 269 clusters detected in the //-band we can determine from a statistical point of 
view both the masses of the young clusters, and what fraction of all the detected infrared 
clusters are young. The main panel of Fig. 6 shows the absolute Mpi6ow magnitude vs. the 
EW of Paa for young clusters selected in the H-band. We show clusters with observed EW 
of Paa > 32 A or ages of less than ~ 7 million years. We also show outputs of Starburst99 
(Leitherer et al. 1999) for clusters with masses M c i uster = 10 4 M and M c i uster = 10 3 M 
formed in an instantaneous burst with a Salpeter IMF between 1 and 100 M . The arrow 
shows what would be the effect of correcting the observed absolute Mpi6ow magnitudes by 
2 magnitudes of visual extinction for the case of a foreground dust screen (assuming extinction 
to stars and line emission is the same). From the colors of bright clusters (Table 1) and the 
color map in Fig. 2 it is clear that the young clusters are not very obscured, and thus the 
observed absolute magnitudes can be used to infer the stellar mass. From Fig. 6 we find that 
the young clusters have stellar masses of between 1O 3 M and 1O 4 M , in good agreement 
with the findings of Harris et al. (2004). Unfortunately for older clusters where we do 
not have EW of Paa based ages, we cannot derive the stellar masses as the near-infrared 
mass-to-light ratio is dependent on the age of the stellar population. 

The insert of Fig. 6 also shows the distribution of absolute M F160 w magnitudes of all 
the clusters detected and that of the the young (age 7 million years) clusters. We find that 
approximately 20 — 30% of the clusters detected in the H-band have ages of less than 7 



4 However, in general for higher metallicity galaxies, Paa must be used in this way with some caution due 
to the delayed episode of ionization from WR stars (Rigby & Rieke 2004). 
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million years. This fraction remains roughly constant in all magnitude bins. 



-30- 



Table 6. Results from the modelling of the SEDs of the nuclear region and galaxy. 



Property 



Nuclear Region 



Entire Galaxy 1 



Star Formation 
Age of stellar population 

Effective optical depth 
Size of the emitting region 
Mass of HI 
Dust Mass 
Stellar Mass 2 
Bolometric Luminosity 



Burst 
t = 3.2 x 10 6 yr 
ry(eff) = 15.7 
r t = 21 pc 
M(H) = 6.2 x 1O 6 M 
M d = 5.0 x 1O 4 M 
M* = 2.8 x 1O 6 M 
L bol = 1.0x 1O 9 L 



Infall model 
t = 30 x 10 6 yr 
ry(eff) = 0.9 
r t = 252 pc 
M(H) = 5.0 x 1O 7 M 
M d = 9.5 x 1O 4 M 
M* = 1.0 x 1O 7 M 
L bol = 1.8xlO 9 L 



Note. — 1 Results from the modelling of the entire galaxy are from Takagi et al. (2003b), 
using an infall model of chemical evolution with a timescale of to — 100 Myr (i.e., SFR 
oc Mg as /to), recomputed for the distance assumed in this paper (d — 4.1 Mpc). 
2 Stellar masses are for a Salpeter IMF with mass cutoffs of M\ ow = 0.1 M and M up = 60 M . 
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Fig. 8. — Spex L-band spectrum of NGC 5253 and comparison galaxies, except NGC 4945 
for which the data are from Spoon et al. (2003) . The spectra have been shifted to rest-frame 
wavelengths. We also mark the most prominent emission features. 
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7. The Nuclear SED of NGC 5253 
7.1. The Model 

We used the SED model of Takagi, Vansevicius & Arimoto (2003) to reproduce the UV 
to mid-infrared SED of the nuclear region of NGC 5253 (Table 4). This model makes use of 
the population synthesis model of Kodama & Arimoto (1997) to predict the stellar emission. 
For the nucleus of this galaxy we assumed an instantaneous burst of star formation with solar 
metallicity and a Salpeter IMF with lower and upper mass cutoffs of M\ ow = 0.1 M and 

— 60 M , respectively. In the SED model of Takagi et al. (2003a), stars are centrally 
concentrated in a spherical region with a density profile given by the King model, while the 
dust grains are homogeneously distributed throughout the emitting region. The equations of 
radiative transfer are numerically solved by considering multiple scattering of photons and 
self-absorption by dust (see Takagi et al. 2003a for more details). The main parameters 
of this model are the age of the stellar population, the equivalent V-band optical depth ry 
measured from the center to the outer edge of the region, and the size of the emitting region, 
r t . Note that ry is defined by the column density of dust, whereas ry (eff) is defined by the 
luminosity ratio of input and output at V-band. The latter optical depth was used to verify 
the SED models of star forming galaxies of Takagi, Arimoto, & Hanami (2003b) because 
it is directly comparable with observations. The extinction curve can be chosen from three 
types: MW, LMC, and SMC. Using this model, we manually found the best-fitting model 
to the nuclear SED, i.e. we calculated a new model in each step of the fitting process until 
convergence was reached. 



7.2. Results from the modelling 

In Table 6 the results of the modelling of the nuclear SED of NGC 5253 are presented 
and compared with those from the modelling of the entire galaxy emission (Takagi et al. 
2003b). The fit to the nuclear SED is shown in Fig. 7. There is a degeneracy in several of 
the model parameters — age, ry, and r t — for the overall nuclear SED from A ~ 1 /im to 
~ 10 /im, but the depth of tbe silicate feature at 9.7 /im can be used to constrain the model 
parameter ranges. We could not obtain reasonable fits with either a LMC or a SMC-type 
extinction curve. For both these extinction curves the fraction of silicate grains is larger 
than that in the MW-type (see Takagi et al. 2003a). We found a conservative limit of 
Ty < 30 from the observed silicate absorption (see also Aitken et al. 1982; Kawara et al. 
1989 and Fig. 7) with the MW-type extinction curve. With this ry, the observed nuclear 
SED from ~ 1 /im to 10 /im was reproduced with SEDs in which the contribution from red 
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supergiants at near-infrared wavelengths is still negligible, i.e. < 4 Myr in the adopted 
population synthesis model. An older age would need values of the optical depth ry > 30 to 
fit the SED from ~ 1 jum to 10 /iin because the relative brightness of the intrinsic SED in the 
near-infrared increases as the system ages. Since the observed depth of the silicate feature 
is relatively shallow, we fit values of the optical depth ry < 20 and the age t < 3.5 Myr. 
The fitted age is in good agreement with the values obtained in Section 5.1. To compare the 
fitted optical depth with observations of the gas emission, we need to use Ty(eff). When ry 
is large, it is not straightforward to compare Ty(eff) with Ay as derived from the emission 
lines. For a meaningful comparison, we need to choose a wavelength to calculate r(eff) 
at which there be still significant stellar emission from the central region and be close to 
the wavelengths used to estimate Ay. For the model with t = 3.2 Myr and Ty = 20, we 
derived an effective optical depth at 2.2 fim of Tx(eff) = 1.4 from the luminosity ratio of the 
input and the attenuated stellar emission. For the MW extinction curve, this corresponds 
to ry(eff) = 15.7 or Ay = 17mag. This is in good agreement with the estimate of Ay from 
the infrared lines in Section 5.2.1. 

Considering these results, the best fit to the nuclear SED is reached for a stellar mass 
of 2.8 x 1O 6 M (for a Salpeter IMF with mass cutoffs M low = 0.1 M and M up = 60 M ) 
formed in an instantaneous burst, with an age of 3.2 Myr, Ty = 20 with the MW extinction 
curve. The fitted size for the emitting region is 21 pc. The stellar mass needed to account 
for the observed nuclear SED of NGC 5253 is in good agreement with the stellar mass of 
the double cluster C1+C2 derived in Section 5.3 (taking into account the correction factor 
to use the same lower mass cutoff of the IMF). The gas and dust masses are 6 x 1O 6 M 
and 5 x 10 4 M , respectively. The resulting bolometric luminosity is 10 9 L . One interesting 
result from comparing the fitting results of the nuclear and galaxy SEDs (Table 6) is that 
up to half of the total dust mass in NGC 5253 appears to be concentrated in the nuclear 
region, with the consequence that a similar fraction of the bolometric luminosity is due to 
the central cluster. 

We find that the fluxes at A < 1 fim are not well reproduced with this model (Fig. 7). 
It is possible that these fluxes are still contaminated by background radiation from the 
underlying galaxy (that is, not arising from the double cluster) and/or photons leaking 
through gas clouds, an effect that is not taken into account in the model. However a more 
likely explanation is the complicated geometry of the central double cluster. As discussed 
in Section 4, most of the UV and optical fluxes measured through the small aperture come 
from cluster CI which is not significantly obscured, at least as infered from its observed near- 
infrared colors (Table 1 and Section 5.2.2) and color maps (Fig. 2). Cluster C2 on the other 
hand, only starts becoming visible in the near-infrared and dominates the nuclear emission 
at wavelengths longer than approximately 2 /iin, indicating that it is deeply embedded. It is 
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clear that the near-infrared to mid-infrared SED is dominated by C2, and thus the UV and 
optical fluxes should be taken as upper limits to the emission of C2. This is suggestive of 
different dust geometry and optical depth for each of the two nuclear clusters. 

The ISO 15 jum flux and the ground-based 20 /im flux observed with the large aperture 
size (Frogel et al. 1982) could not be well reproduced with the model, indicating that there 
may be some contribution from emitting sources other than the central double cluster. Note 
however, that the smaller aperture fluxes of Lebofsky & Rieke (1979) and Gorjian et al. 
(2001) are well reproduced with the model. At A < 20 /im, the emission from very small 
grains (graphite) dominates, and because of this, there is a 'shoulder' in the SED. Although 
the small aperture 20 /im data seems to indicate a single peak of hot dust emission, no SED 
model has such a single peak. If we reduced the size of the emitting region to move the peak 
of the SED toward shorter wavelengths, then the near-infrared fluxes would be significantly 
over-estimated. Therefore, there should be a true peak at 30 /im due to large grains, unless 
the emission from very small grains is significantly deficient. Finally, the best fit model 
shows that the 3.3 /im polycyclic aromatic hydrocarbon (PAH) feature is still present, but 
is not very prominent since the continuum produced by very small grains is high, effectively 
reducing the equivalent width of this feature. A similar fit to the nuclear SED of NGC 5253 
could be obtained if we reduced the amount of PAH carrier, since the infrared SED remains 
the same, except for the PAH feature (i.e., PAHs have negligible contribution to continuum). 
We return to this issue in Section 7.3. 



7.3. The non-detection of the 3.3 /im PAH feature in the nucleus of 

NGC 5253: age or metallicity effect? 

Most massive star-forming galaxies display infrared emission features at 3.3, 6.2, 7.7, 
8.6, and 11.3 /im (Roche et al. 1991; Uchida et al. 2000; Helou et al. 2000) that have been 
attributed to bending and stretching modes of CH and CC bonds in PAHs (Duley & Williams 
1981; Leger & Puget 1984; Puget & Leger 1989). Among others, Roche et al. (1991) and 
Clavel et al. (2000) have shown that the strength of the 11.3 /im PAH feature depends 
strongly on environment: galaxies classified as Seyfert 1 show little to no 11.3 /im emission, 
while many type 2 AGN (see also Imanishi, 2002 and 2003 for a study of the 3.3 /im feature) 
and most pure starburst galaxies show strong PAH features. This behavior implies that as 
the radiation field density increases, the PAH carrier is progressively destroyed. This picture 
is supported by M82 where the 3.3 /im PAH emission is found to be anti-correlated with 
both the centrally peaked continuum emission and also with the Br7 emission (Normand et 
al. 1995). A similar situation is observed in NGC 253 where there is a decrease of the ratio 
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of 3.3 /im PAH emission to continuum at the sites of strong star formation (Tacconi-Garman 
et al. 2004). 

Dwarf galaxies in the metallicity range 12 + logO/H = 7.6 — 8.6 tend to show weak or no 
mid-infrared PAH features from ISO data (Madden 2002; Thuan, Sauvage, & Madden 1999) 
or ground-based observations (Roche et al. 1991), including NGC 5253 with 12 + logO/H = 
8.4. We are obtaining L-band spectroscopy of a sample of starburst galaxies with a range 
of metallicities using Spex on the IRTF. Fig. 8 are the spectra for some of these galaxies, 
including a second dwarf galaxy, He2-10, with an age of ~ 5 — 6 Myr and metallicity of 
12 + logO/H = 8.9 (Kobulnicky, Kennicutt, & Pizagno 1999). For NGC 2903 we have 
spectra for the quiescent nucleus and one of the bright Hll regions to the north of the 
nucleus (see Alonso-Herrero et al. 2001). The 3.3 /xm feature is clearly detected in He2-10 
and in the bright Hll region of NGC 2903, as well as in all starburst galaxies in our sample 
with metallicities 12 + log O/H > 8.9. We also show that the highly obscured starburst/AGN 
NGC 4945 (L-band spectroscopy from Spoon et al. 2003) has bright 3.3 /im PAH emission, 
indicating that this feature can be detected regardless of the extinction (Imanishi 2002). 

The effects of metallicity and a harsh environment are coupled for young starbursts 
as in the nucleus of NGC 5253. As metallicity decreases stars present harder spectra, and 
Wolf-Rayet stars require larger mass progenitors thus creating a hard and intense radiation 
field that could suppress the formation of the 3.3 /im PAH feature. As discussed in Tacconi- 
Garman et al. (2004) and references therein, in strong radiation fields the efficiency of 3.3 
PAH emission can be lowered by photoionization of the PAH molecules, or by destruction 
of the PAH molecules. Rigby & Rieke (2004) have demonstrated that the observed mid- 
infrared fine structure lines of the low metallicity starbursts NGC 5253 and II Zw40 can only 
be reproduced with a Salpeter IMF extending to 40 — 60 M , whereas for their sample of high 
mass, solar metallicity starbursts the interstellar radiation field will not include the outputs 
of stars above 40 M Q because most of them are probably embedded in compact H II regions. 
Although the number of low metallicity starbursts analyzed in Rigby & Rieke (2004) is small, 
there seems to be an indication that in these systems the interstellar medium is less effective 
in confining ultracompact Hll regions, and thus a harsher radiation field environment is 
present that can effectively decrease or even suppress the 3.3 /im PAH emission. 



8. Discussion and Conclusions 

We have analyzed iLST/NICMOS observations (1.1 — 2.2 /tm) and obtained 1.9 — 4.1 /tm 
spectroscopy of the central region of the dwarf galaxy NGC 5253. The NICMOS images have 
revealed the presence of a double cluster (C1+C2) in the nucleus of the galaxy separated by 
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0.3 — 0.4" or 6 — 8pc for the assumed distance of d — 4.1 Mpc. This double cluster is also 
a bright double source of Pact emission. We have also analyzed HST/WFPC2 observations 
and found that the western cluster (C2) is almost entirely obscured at UV and optical 
wavelengths, but becomes the brightest source in the galaxy at A > 2 /zm. The double 
cluster C1+C2 appears to be coincident with the double radio nebula detected at 1.3 cm by 
Turner et al. (2000), and it is likely to be responsible for most of the mid-infrared nuclear 
emission in this galaxy (see also Turner & Beck 2004). 

The high spatial resolution of the HST/N1CMOS images allows us to measure the near- 
infrared magnitudes, EW of Pact and Pact luminosities of CI and C2. If the star formation 
in the double cluster C1+C2 occurred in an instantaneous burst, CI and C2 have virtually 
the same age (3 — 4 million years). This young age for the double cluster is consistent with 
estimates by other authors for cluster CI from UV/optical imaging and spectroscopy (e.g., 
Calzetti et al. 1997; Harris et al. 2004), and the presence of Wolf-Rayet features in the 
optical spectrum of the nucleus of the galaxy (Walsh & Roy 1987; Conti 1991). Although 
both nuclear clusters appear to be coeval, C2 is more massive and obscured than CI. The 
stellar mass for C2 is in the range of M C2 = 7.7 x 10 5 - 2.6 x 10 6 M (for a Salpeter IMF in 
the mass range 0.1 — 100 M Q ), depending on the dust emission contribution to the observed 
X-band magnitude, thus putting C2 in the category of super star clusters. The stellar mass 
of CI is Mci ^ 1.3 x 1O 5 M . 

The fact that both CI and C2 appear to be coeval but C2 is much more massive than 
CI appears to conflict with the naive expectation that more massive clusters would blow 
away their birth clouds first. Using K and L-band hydrogen recombination lines and near- 
infrared colors we have inferred a very high obscuration for C2, suggesting that this young 
star cluster may be still (partially) embedded in its natal cloud. The age derived using the 
EW of Pact assumes that the extinction to the stars and gas is the same. If this is not the 
discussed in Calzetti et al. (1997), that is, stars, gas and dust are mixed, then the 
derived age of C2 would be only an upper limit. 

Turner & Beck (2004) have put forward some alternative explanations to that of the 
presence of a coeval double cluster in the nucleus of NGC 5253. One possibility they discuss 
is that cluster CI is not a cluster but a reflection nebula from a gap opening in the cocoon 
around the infrared cluster. Their second scenario is that the supernebula is cometary and 
that what we are seeing in the nucleus of NGC 5253 is motion of the H II region/ cluster to the 
west and north. We favor the interpretation of a double cluster in the nucleus of NGC 5253 
that has resulted from triggered star formation, due to infalling gas along the minor axis 
of the galaxy (see Meier et al. 2002). The south-east elongation of the radio supernebula 
detected by Turner & Beck (2004) shows a morphology similar to the Pact emission (Fig. 4, 
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left panel), but we find that the bright nuclear Pact emission is extended over scales of a 
few arcseconds, whereas the radio supernebula is very compact (~ 0.05", see Turner & Beck 
for more details). From the if-band NICMOS image of the central 3" x 3" it is clear that 
there are a number of near-infrared clusters with no Pact emission, but there must be other 
ionizing clusters (not detected as bright sources of 1.6 /zm emission) that are responsible for 
the extended Pact emission. Moreover, we have derived a very young age (3 — 4 million 
years) using the i^-band spectrum for the circumnuclear region excluding CI and C2. This 
all seems to indicate that what we are seeing are the effects of self-propagating star formation 
in the central region of NGC 5253. 

In addition to the nuclear double cluster, we have identified a total of 269 star clusters 
in the if-band over the observed region of 19" x 19", and measured their near-infrared 
magnitudes as well as the EW of Pact. Assuming that the EW of hydrogen recombination 
lines are good indicators of the age, we find that 20 — 30% of the detected clusters in the 
if -band have ages younger than approximately 7 million years. These young clusters have 
stellar masses of between 3 x 10 3 M© and 3 x 1O 4 M . 

We finally model the nuclear UV to mid-infrared SED of NGC 5253 using the model 
of Takagi et al. (2003a), and compare it to the SED of the entire galaxy. For the nuclear 
SED we have taken special care to measure and subtract the underlying galaxy emission to 
isolate as much as possible the UV to mid-infrared emission of the double star cluster. For 
a Salpeter IMF we find that the nuclear infrared SED is well reproduced with a total stellar 
mass of M* = 2.8 x 1O 6 M and an age of 3.2 million years in good agreement with the 
estimates using the photometry and spectroscopy of C1+C2. The modelling of the nuclear 
SED shows that the young starburst is very obscured, Ay = 17 mag, in good agreement 
with extinction estimates obtained from near-infrared hydrogen recombination lines. The 
UV-optical SED of the nuclear region is not well reproduced by the model probably due to 
the complicated dust and emitting source geometry of the nuclear region. The comparison 
between the nuclear and entire galaxy SED modelling shows that up to 50% of the dust mass 
in this galaxy is concentrated in the nucleus. 

The model of the nuclear SED of NGC 5253 also predicts a moderately bright 3.3 fxm 
PAH feature that is not observed in our nuclear L-band spectrum. The 3.3 /iin PAH feature 
along with the mid-infrared PAH features are detected in most massive starburst galaxies, 
but appear to be weak or absent in low mass and metallicity dwarf star forming galaxies. 
Rigby & Rieke (2004) have recently demonstrated that the star formation properties of 
NGC 5253, in particular the mid-infrared fine structure line ratios, can only be reproduced 
with a Salpeter IMF where the upper mass cutoff is M up > 40 — 60 M Q . In the other more 
massive and metal rich starbursts in their sample the interstellar radiation field may not 
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include the outputs of stars above 40 M because most of them are probably embedded in 
compact Hll regions. The fact that low metallicity stars produce a harder radiation field 
together with the presense of massive stars in the nuclear starburst of NGC 5253 result in a 
more intense and harder radiation field that could effectively suppress the production of the 
3.3 /iin PAH feature. 

We thank Henrik Spoon for providing us with the L-band spectrum of NGC 4945, and 
J. Turner and S. Beck for providing us with the VLA radio data of NGC 5253. We are also 
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NASA Contract 960785 through the Jet Propulsion Laboratory, and the Spanish Programa 
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